A role for magnetic susceptibility in synthetic computed tomography.

PURPOSE: Radiotherapy treatment planning based on magnetic resonance imaging (MRI) benefits from increased soft-tissue contrast and functional imaging. MRI-only planning is attractive but limited by the lack of electron density information required for dose calculation, and the difficulty to differentiate air and bone. MRI can map magnetic susceptibility to separate bone from air. A method is introduced to produce synthetic CT (sCT) through automatic voxel-wise assignment of CT numbers from an MRI dataset processed that includes magnetic susceptibility mapping. METHODS: Volumetric multi-echo gradient echo datasets were acquired in the heads of five healthy volunteers and fourteen patients with cancer using a 3 T MRI system. An algorithm for CT synthesis was designed using the volunteer data, based on fuzzy c-means clustering and adaptive thresholding of the MR data (magnitude, fat, water, and magnetic susceptibility). Susceptibility mapping was performed using a modified version of the iterative phase replacement algorithm. On patient data, the algorithm was assessed by direct comparison to X-ray computed tomography (CT) scans. RESULTS: The skull, spine, teeth, and major sinuses were clearly distinguished in all sCT, from healthy volunteers and patients. The mean absolute CT number error between X-ray CT and sCT in patients ranged from 78 and 134 HU. CONCLUSION: Susceptibility mapping using MRI can differentiate air and bone for CT synthesis. The proposed method is automated, fast, and based on a commercially available MRI pulse sequence. The method avoids registration errors and does not rely on a priori information, making it suitable for nonstandard anatomy.

Microdosimetric Evaluation of Current and Alternative Brachytherapy Sources-A Geant4-DNA Simulation Study.

PURPOSE: Radioisotopes such as 75Se, 169Yb, and 153Gd have photon energy spectra and half-lives that make them excellent candidates as alternatives to 192Ir for high-dose-rate brachytherapy. The aim of the present study was to evaluate the relative biological effectiveness (RBE) of current (192Ir, 125I, 103Pd) and alternative (75Se, 169Yb, 153Gd) brachytherapy radionuclides using Monte Carlo simulations of lineal energy distributions. METHODS AND MATERIALS: Brachytherapy sources (microSelectron v2 [192Ir, 75Se, 169Yb, 153Gd], SelectSeed [125I], and TheraSeed [103Pd]) were placed in the center of a spherical water phantom with a radius of 40 cm using the Geant4 Monte Carlo simulation toolkit. The kinetic energy of all primary, scattered, and fluorescence photons interacting in a scoring volume were tallied at various depths from the source. Electron tracks were generated by sampling the photon interaction spectrum and tracking all the interactions down to 10 eV using the event-by-event capabilities of the Geant4-DNA models. The dose mean lineal energy (y¯D) values were obtained through random sampling of transfer points and overlaying spherical scoring volumes within the associated volume of the tracks. The scoring volume diameter was determined by fitting the y¯D ratio for 125I to its observed RBE. RESULTS: y¯D increased with the increasing distance from the source for 192Ir, 75Se, and 169Yb, remained constant for 153Gd and 125I, and decreased for 103Pd. The diameter at which the y¯D ratio coincided with the RBE of 1.15 to 1.20 for 125I was ∼25 to 40 nm. The RBE (reference 1 MeV photons) at high doses and dose rates for 192Ir, 75Se, 169Yb, 153Gd, 125I, and 103Pd was 1.028 to 1.034, 1.05 to 1.07, 1.12 to 1.15, 1.16 to 1.21, 1.15 to 1.20, and 1.17 to 1.22, respectively. CONCLUSIONS: The radiation quality of the radionuclides under investigation was greater than that of high-energy photons. The present study has provided a set of values to modify the prescription doses for brachytherapy to account for the variation in radiation quality among radionuclides.

Microdosimetry calculations for monoenergetic electrons using Geant4-DNA combined with a weighted track sampling algorithm.

The aim of this study was to calculate microdosimetric distributions for low energy electrons simulated using the Monte Carlo track structure code Geant4-DNA. Tracks for monoenergetic electrons with kinetic energies ranging from 100 eV to 1 MeV were simulated in an infinite spherical water phantom using the Geant4-DNA extension included in Geant4 toolkit version 10.2 (patch 02). The microdosimetric distributions were obtained through random sampling of transfer points and overlaying scoring volumes within the associated volume of the tracks. Relative frequency distributions of energy deposition f(>E)/f(>0) and dose mean lineal energy ([Formula: see text]) values were calculated in nanometer-sized spherical and cylindrical targets. The effects of scoring volume and scoring techniques were examined. The results were compared with published data generated using MOCA8B and KURBUC. Geant4-DNA produces a lower frequency of higher energy deposits than MOCA8B. The [Formula: see text] values calculated with Geant4-DNA are smaller than those calculated using MOCA8B and KURBUC. The differences are mainly due to the lower ionization and excitation cross sections of Geant4-DNA for low energy electrons. To a lesser extent, discrepancies can also be attributed to the implementation in this study of a new and fast scoring technique that differs from that used in previous studies. For the same mean chord length ([Formula: see text]), the [Formula: see text] calculated in cylindrical volumes are larger than those calculated in spherical volumes. The discrepancies due to cross sections and scoring geometries increase with decreasing scoring site dimensions. A new set of [Formula: see text] values has been presented for monoenergetic electrons using a fast track sampling algorithm and the most recent physics models implemented in Geant4-DNA. This dataset can be combined with primary electron spectra to predict the radiation quality of photon and electron beams.

Proton and light ion RBE for the induction of direct DNA double strand breaks.

PURPOSE: To present and characterize a Monte Carlo (MC) tool for the simulation of the relative biological effectiveness for the induction of direct DNA double strand breaks (RBEDSB (direct)) for protons and light ions. METHODS: The MC tool uses a pregenerated event-by-event tracks library of protons and light ions that are overlaid on a cell nucleus model. The cell nucleus model is a cylindrical arrangement of nucleosome structures consisting of 198 DNA base pairs. An algorithm relying on k-dimensional trees and cylindrical symmetries is used to search coincidences of energy deposition sites with volumes corresponding to the sugar-phosphate backbone of the DNA molecule. Strand breaks (SBs) are scored when energy higher than a threshold is reached in these volumes. Based on the number of affected strands, they are categorized into either single strand break (SSB) or double strand break (DSB) lesions. The number of SBs composing each lesion (i.e., its size) is also recorded. RBEDSB (direct) is obtained by taking the ratio of DSB yields of a given radiation field to a (60)Co field. The MC tool was used to obtain SSB yields, DSB yields, and RBEDSB (direct) as a function of linear energy transfer (LET) for protons ((1)H(+)), (4)He(2+), (7)Li(3+), and (12)C(6+) ions. RESULTS: For protons, the SSB yields decreased and the DSB yields increased with LET. At ≈24.5 keV µm(-1), protons generated 15% more DSBs than (12)C(6+) ions. The RBEDSB (direct) varied between 1.24 and 1.77 for proton fields between 8.5 and 30.2 keV µm(-1), and it was higher for iso-LET ions with lowest atomic number. The SSB and DSB lesion sizes showed significant differences for all radiation fields. Generally, the yields of SSB lesions of sizes ≥2 and the yields of DSB lesions of sizes ≥3 increased with LET and increased for iso-LET ions of lower atomic number. On the other hand, the ratios of SSB to DSB lesions of sizes 2-4 did not show variability with LET nor projectile atomic number, suggesting that these metrics are independent of the radiation quality. Finally, a variance of up to 8% in the DSB yields was observed as a function of the particle incidence angle on the cell nucleus. This simulation effect is due to the preferential alignment of ion tracks with the DNA nucleosomes at specific angles. CONCLUSIONS: The MC tool can predict SSB and DSB yields for light ions of various LET and estimate RBEDSB (direct). In addition, it can calculate the frequencies of different DNA lesion sizes, which is of interest in the context of biologically relevant absolute dosimetry of particle beams.

On the consistency of Monte Carlo track structure DNA damage simulations.

PURPOSE: Monte Carlo track structures (MCTS) simulations have been recognized as useful tools for radiobiological modeling. However, the authors noticed several issues regarding the consistency of reported data. Therefore, in this work, they analyze the impact of various user defined parameters on simulated direct DNA damage yields. In addition, they draw attention to discrepancies in published literature in DNA strand break (SB) yields and selected methodologies. METHODS: The MCTS code Geant4-DNA was used to compare radial dose profiles in a nanometer-scale region of interest (ROI) for photon sources of varying sizes and energies. Then, electron tracks of 0.28 keV-220 keV were superimposed on a geometric DNA model composed of 2.7 × 10(6) nucleosomes, and SBs were simulated according to four definitions based on energy deposits or energy transfers in DNA strand targets compared to a threshold energy ETH. The SB frequencies and complexities in nucleosomes as a function of incident electron energies were obtained. SBs were classified into higher order clusters such as single and double strand breaks (SSBs and DSBs) based on inter-SB distances and on the number of affected strands. RESULTS: Comparisons of different nonuniform dose distributions lacking charged particle equilibrium may lead to erroneous conclusions regarding the effect of energy on relative biological effectiveness. The energy transfer-based SB definitions give similar SB yields as the one based on energy deposit when ETH ≈ 10.79 eV, but deviate significantly for higher ETH values. Between 30 and 40 nucleosomes/Gy show at least one SB in the ROI. The number of nucleosomes that present a complex damage pattern of more than 2 SBs and the degree of complexity of the damage in these nucleosomes diminish as the incident electron energy increases. DNA damage classification into SSB and DSB is highly dependent on the definitions of these higher order structures and their implementations. The authors' show that, for the four studied models, different yields are expected by up to 54% for SSBs and by up to 32% for DSBs, as a function of the incident electrons energy and of the models being compared. CONCLUSIONS: MCTS simulations allow to compare direct DNA damage types and complexities induced by ionizing radiation. However, simulation results depend to a large degree on user-defined parameters, definitions, and algorithms such as: DNA model, dose distribution, SB definition, and the DNA damage clustering algorithm. These interdependencies should be well controlled during the simulations and explicitly reported when comparing results to experiments or calculations.

Monte Carlo role in radiobiological modelling of radiotherapy outcomes.

Radiobiological models are essential components of modern radiotherapy. They are increasingly applied to optimize and evaluate the quality of different treatment planning modalities. They are frequently used in designing new radiotherapy clinical trials by estimating the expected therapeutic ratio of new protocols. In radiobiology, the therapeutic ratio is estimated from the expected gain in tumour control probability (TCP) to the risk of normal tissue complication probability (NTCP). However, estimates of TCP/NTCP are currently based on the deterministic and simplistic linear-quadratic formalism with limited prediction power when applied prospectively. Given the complex and stochastic nature of the physical, chemical and biological interactions associated with spatial and temporal radiation induced effects in living tissues, it is conjectured that methods based on Monte Carlo (MC) analysis may provide better estimates of TCP/NTCP for radiotherapy treatment planning and trial design. Indeed, over the past few decades, methods based on MC have demonstrated superior performance for accurate simulation of radiation transport, tumour growth and particle track structures; however, successful application of modelling radiobiological response and outcomes in radiotherapy is still hampered with several challenges. In this review, we provide an overview of some of the main techniques used in radiobiological modelling for radiotherapy, with focus on the MC role as a promising computational vehicle. We highlight the current challenges, issues and future potentials of the MC approach towards a comprehensive systems-based framework in radiobiological modelling for radiotherapy.